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Universal relations in the finite-size correction terms of two-dimensional Ising models 
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Quite recently, Izmailian and Hu [Phys. Rev. Lett. 86, 5160 (2001)] studied the finite-size 
correction terms for the free energy per spin and the inverse correlation length of the critical two- 
dimensional Ising model. They obtained the universal amplitude ratio for the coefficients of two 
series. In this study we give a simple derivation of this universal relation; we do not use an explicit 
form of series expansion. Moreover, we show that the Izmailian and Hu's relation is reduced to a 
simple and exact relation between the free energy and the correlation length. This equation holds 
at any temperature and has the same form as the finite-size scaling. 
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Universality and scaling are two basic concepts in the 
study of phase transitions and critical phenomena [|l],p| . 
The critical properties are universal to a large extent 
depending on a few parameters, such as the space di- 
mensionality and the symmetry of the order parameter. 
Critical exponents, critical amplitude ratios, and scal- 
ing functions are examples of universal quantities P,H. 
Finite-size scaling [^|y] has been increasingly important, 
partly due to the progress in the theoretical understand- 
ing of finite-size effects, and partly due to the application 
to the analysis of simulational results. Recently, more 
attention has been paid to the universality of finite-size 
scaling functions M for both percolation models S and 
Ising models Q . In two dimensions the relevance of the 
finite-size properties to the conformal field theory is an- 
other source of interest jlO| . 

Quite recently, Izmailian and Hu 111] studied the finite- 
size correction terms for the free energy per spin and 
the inverse correlation length of critical two-dimensional 
(2D) Ising models fl^-p^. They obtained the universal 
amplitude ratio for the coefficients of two series. Let us 
denote the free energy per spin and the inverse correla- 
tion length for A^ X cxD lattice as /n and ^^^ , respectively. 
Then, Izmailian and Hu pd| obtained analytic expres- 
sions for the finite-size correction coefficients Uk and bk 
defined by 
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at the critical! ty for the square (sq), honeycomb (he), 
and plane-triangular (pt) lattices. Here, foo denotes the 
value of /jv as A^ — > oo. They found that 
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for all of these lattices, that is, the amplitude ratio bk/ak 
is universal. They also obtained similar expansions for 



the critical ground state energy E^°^ and the first energy 
gap {E^^'' — i?(°)) of the one-dimensional (ID) quantum 
XY model with uniform field ||T^ at the critical field, and 
found that the amplitude ratios of two coefficients have 
the same values as Eq. (^). 

The lowest-order correction terms ai and 6i are re- 
lated to the central charge c and the magnetic scaling 
field xh in the conformal field theory by ai = ctt/6 
and 6i — 2ttxh, respectively. The finite-size correction 
terms of 2D systems have also received current interest. 
Quieroz |16| numerically studied the finite-size correc- 
tions to scaling of correlation lengths and free energies 
of the critical 2D Ising and 3-state Potts models. The 
finite-size corrections to the energy and specific heat of 
the critical 2D Ising model were also analyzed by Salas 
pTf . Moreover, Salas and Sokal studied several universal 
amplitude ratios for the critical 2D Ising model |18[ |. 

In this paper we will give a simple derivation of the 
universal amplitude ratio, Eq. (g). Instead of using an 
explicit form of expansion in 1/A'^, which is available only 
at the criticality, we show a finite-size-scaling like relation 
which holds at an arbitrary temperature. The universal 
amplitude ratio can be readily derived from this relation. 

Let us start with comparing the free energy per spin 
and the inverse correlation length for size N and N/2. 
Assuming the expansion, Eq. dl^), we have the expression 
for the difference of Jn and /jv/2 as 
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In a similar way, using Eq. (g), we have 
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From Eqs. (|j) and (|^), we get the following statement: 
If the relation 
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is satisfied, we have the relation of two coefficients, 
Eq. (^. If we can show Eq. ^ directly, explicit ex- 
pressions for the coefficients Ofe and bk are not necessary. 
Equation (||) is compatible to the finite-size scaling [g|J|] 
of the singular part of the free energy for N x oo sys- 
tems, /n oc {1/N)S,]^^ {N -^ cxj); namely, Eq. (||) should 
hold asymptotically at the critical point. In fact, as we 
show in the following, Eq. (|6|) holds exactly at any tem- 
perature for several models that belong to the 2D Ising 
universality class. 

First, we deal with the 2D Ising model defined by the 
Hamiltonian 
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where the Ising variable Si takes ±1, f3 = {kBT)~^, and 
the summation is taken over the nearest-neighbor pairs of 
sites for N x oo lattices. Using a transfer matrix method 
|^,|0[, we can calculate the free energy per spin, /jv, and 
the inverse correlation length, ^^^ , through the relations 



/Ar = — InAo, 
G' = ^ln(Ao/Ai), 



(8) 
(9) 



where Aq and Ai are the largest and the second largest 
eigenvalues of the transfer matrix. Here, a geometric fac- 
tor ( is 1, 2/\/3 and l/-\/3 for sq, he, and pt lattices, 
respectively. Exact expressions for eigenvalues Aq and 
Ai are available for sq |l|,|l|,|§l, he Q, and pt |l|] 
lattices. 

Let us consider the Onsager solution for the sq lattice 
with periodic boundary conditions |lj] . The two leading 
eigenvalues of the transfer matrix are given by 
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Then, we get the free energy per spin f^, Eq. (||), and 
the inverse correlation length ^^^ , Eq. (^) , as 

1 1 ^"^ 

/at = - ln(2 sinh 2 J) + — ^ 72^+1, (13) 
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To express fN/2 and ^n/2, we replace N by N/2 and 
Ir by 7^ in Eqs. (O) and (Hj), respectively. Here 7^ is 



the value for size N/2, and is related to 7^ through the 
relation 



, , cosh^ 2 J rn 

cosh 7^ = . ^ ^ ^ — cos -^^^-j^ ~ cosh 72^ 



sinh 2 J 



N/2 



that is. 



72r 



(15) 



(16) 



Equation (16) is the basic relation for our argument. Us- 
ing the relation (|l(^), we have 

fN/2 ~ In = ^{2(72 + 76 H y- I2N-2) 

-(71+73H l-72jv-i)|- (17) 

Similarly, we have 
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Comparing Eqs. (|lj) and (|18|), we arrive at the desired 
relation, Eq. (|g). This implies that the ratio for the finite- 
size correction coefficients bk/ak is given by Eq. (0). 

Next we consider the Ising model on the he lattice. 
The two leading eigenvalues of the transfer matrix are 
given by @ 

Af/2-l 

Ao = (2sinh2J)^exp( ^ 72.+1), (19) 
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In this case, 7^ is implicitly given by 

o nr 
cosh7r — cosh 2J cosh 2 J* — sin — 

N 

1 /o 
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where J* is defined by (cosh2J— l)(cosh2J*-l) — 1. For 
size N/2, we replace N by N/2 and 7^. by 7^ in Eqs. ( pj| ) 
and (I20I). The relation of 7^, the value for size N/2, and 
7r can be obtained by using Eq. (p^); we have the same 
relation as Eq. (|lj), 7^ — 72.- Starting from Eqs. (|^) and 
(P) together with Eqs. ( p^ ) and (pO[), with some algebra, 
we get 
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Thus, we again obtain the relation (ra). In this way, we 
have shown that the aniphtude ratio hk/au is given by 
Eq. (^). For the pt lattice, we may use the star-triangle 
transformation |^1|] to show Eq. (|g) from the result of the 
he lattice without making an explicit calculation. 

Some comments should be added here. We have shown 
the universality of the amplitude ratio, Eq. (0), without 
using the condition of the criticality. This means that 
the property of this universality holds not only at the 
critical point but also at any temperature in 2D Ising 
models. The ratio bk/ak takes the universal value for 
all the temperatures. However, we should note that the 
expression for the inverse correlation length, Eq. (||), is 
valid only for T > T^. For T < Tc, Aq and Ai are de- 
generate in the thermodynamic limit, N —> 00, which 
means the existence of the long-range order M,pO[ . We 
should subtract the long-range order contribution from 
the correlation function when considering the correlation 
length. We may regard bk as the correction amplitude 
for the right-hand side of Eq. (|) for T < T^. 

Izmailian and Hu ||l^ also studied another model 
which belongs to the 2D Ising universality class. The 
ID quantum XY model with uniform field, whose Hamil- 
tonian is given by 
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was exactly solved by Katsura ||l5|]. Here, u^ , a^ and 
tj^ are the Pauli matrices. For < 7 < 1, there is a 
critical magnetic field he = 1, and the phase transition 
of this model belongs to the 2D Ising universality class. 
For 7 = 1 it is also called the ID transverse Ising model. 
The Hamiltonian of Eq. (23) is diagonalized by a Jordan- 



Wigner transformation as 
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where ry^., r]k are fermionic creation and annihilation op- 
erators and 



A(fc) = ^(cos k + h)"^ + (7 sin ky 



(25) 



We should note that the choice of k depends on the 
boundary condition. The ground state energy E^*^' and 
the first energy gap AE — E^^^i — _E(°) of the quantum 
spin model, respectively, correspond to the free energy 
and inverse correlation length for the Ising model, that 
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For the N quantum spin systems with the periodic 
boundary condition, we have [p2|,p3[ 
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where we have used the notation 7^ = K{rTi/N). Let us 
consider the energy for size N and that for N/2 as in 
the classical 2D Ising model. For size N/2, we replace N 
by N/2 and 7^ by 7; in Eqs. (||) and ^. Here, 7^ is 
again the value for N/2, and the relation of 7^ and 7^ 
can be obtained by using Eq. (|2^). As a result, we have 
the same relation as Eq. (K^), that is. 
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Then, we obtain the expression for the difference between 
i?l?) and E^!!l as 
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For the first energy gap, we can also calculate the differ- 
ence between the values for N and N/2. Using Eq. ([2^), 
with some algebra, we finally have 
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From Eqs. (|29|) and (30), we have 



(-^Af/2 ~ ^N ) ^ ^En/2 — AEn, 



(31) 
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which shows that the amplitude ratio of the finite-size 
correction coefficient for the ground state energy (— _E(°)), 
ak, and that for the first energy gap AE, bk, is given by 
Eq. (|). 

The amplitude ratio bk/ak takes the same value as the 
classical case. And this is not only at the critical field 
he — 1 but also for all the magnetic fields. However, it 
should be mentioned that for low field, h < he, Ej^ does 
not necessarily give the first excited energy level. We 
should interpret that bk is the correction amplitude for 
the right-hand side of Eq. [27) ioi h < he- 

To summarize, we have given a simple derivation of the 
universal amplitude ratio, Eq. (0). We have shown the 
relation of Eq. m) and its quantum counterpart, Eq. (pn|). 
Although these equations have the same form as the 
finite-size scaling [§,§, they are exact and valid for an 



arbitrary temperature. One can perform an analytic cal- 
culation of each correction amplitude only at the criti- 
cality [Q. Our key relations are Eqs. (16) and ([2^), and 
we have not used the condition of the critical point. 

Izmailain and Hu pl| used a perturbed conformal field 
theory to understand the correction terms. It is interest- 
ing to study directly the amplitude ratio by using the 
conformal field theory. The present study may give a 
hint for such a direction. 

The finite-size properties depend on the boundary con- 
ditions. Recently, the effect of peculiar boundary condi- 
tions, such as the Mobius strip and the Klein bottle, was 
studied for the 2D Ising model [|4|J2|]. The effect of 
boundary conditions on the finite-size correction terms is 
an interesting subject to study, which is now in progress. 

In this paper, we have started from the expansion of 
the free energy and the inverse correlation length in odd 
powers of 1/A^ in Eqs. (|l|) and (|2|). Away from the crit- 
ical point, correction terms in even powers of 1/A^ may 
appear. Our argument is easily extended to such a case, 
and the main conclusion remains the same. 
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